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We have studied the electronic structure of InN and GaN employing Go Wo calculations based on 
exact-exchange density-functional theory. For InN our approach predicts a gap of 0.7 eV. Taking the 
Burnstein-Moss effect into account, the increase of the apparent quasiparticle gap with increasing 
electron concentration is in good agreement with the observed blue shift of the experimental optical 
absorption edge. Moreover, the concentration dependence of the effective mass, which results from 
the non-parabolicity of the conduction band, agrees well with recent experimental findings. Based 
on the quasiparticle band structure the parameter set for a 4 x 4 kp Hamiltonian has been derived. 
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The group Ill-nitrides A1N, GaN and InN and their 
alloys have become an important class of semiconductor 
materials, in particular for use in optoelectronic devices 
such as green and blue light emitting diodes (LEDs) and 
lasers. Among the three materials InN is still the least 
explored, due to difficulties in synthesizing high quality 
single crystals. Only very recently these problems have 
been overcome .1] , but many of the key band parameters 
have not been conclusively determined until now Q, • 
The most controversially discussed parameter is currently 
still the fundamental band gap of InN. For many years it 
was believed to be approximately 1.9 eV, but essentially 
any value between 0.65 and 2.3 eV has been reported 
in the literature over the last 30 years pj. However, 
more recent experiments on high quality samples grown 
by molecular beam epitaxy (MBE) and recent ab initio 
calculations support a significantly lower value around 
0-7 evil! 111!. 

Different hypotheses have been proposed to explain the 
large variation in the measured band gaps. Defects could 
be responsible for inducing states in the band gap or give 
rise to a pronounced Burnstein-Moss effect due to a shift 
in the Fermi level caused by a high intrinsic electron den- 
sity. Non stoichiometry may increase the defect concen- 
tration or alter the crystaline structure. The formation 
of oxides and oxynitrides would increase the band gap, 
whereas the precipitation of In clusters leads to addi- 
tional features in optical absorption spectra Q, @- In 
this article we demonstrate that first principles calcula- 
tions can contribute to the solution of this fundamental 
question. By combining density-functional theory (DFT) 
with many-body perturbation theory in the GqWq ap- 
proximation , which is currently the method of choice 
for calculating quasiparticle excitations in solids 0, 0] , 
we combine atomistic control over the material with ac- 
curate calculations for the band structure and the band 
gap of stoichiometric and defect free structures. 

Previous ab initio studies were aggravated by the fact 
that DFT calculations in the local-density approximation 
(LDA) predict InN to be metallic in the zinc blende and 
wurtzite structures. Subsequent GqWq calculations only 
open the gap to 0.02 - 0.05 eV [ll [13, while adding 
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FIG. 1: Conduction band of wurtzite InN and GaN aligned 
at the bottom of the conduction band: the circles are the 
Go Wo results, the solid lines the k-p fit using Eq. Q and 
the dashed lines the effective mass band. The inset shows the 
band structure of wurtzite InN. 



self-interaction corrections to the DFT calculations, ei- 
ther in the screened-exchange 0, exact-exchange opti- 
mized effective potential approach (OEPx) wi. or self- 
interaction corrected (SIC) LDA approach [la, 01 yield 
a semiconductor with a band gap of 0.8, 1.0 and 1.6 
eV for the wurtzite phase, respectively. Here we ap- 
ply the Go Wo corrections to OEPx ground state calcu- 
lations, which are fully self-interaction free. We have 
previously shown that this approach yields band gaps for 
II- VI compounds and GaN in very good agreement with 
experiment pl| and for wurtzite InN our value of 0.7 eV 
(see Tab. |lj strongly supports the recent experimental 
findings 0- LJ 13 ■ Similar conclusions were drawn from 
previous Go Wo calculations applied to SIC-LDA ground 
states, however, only after adjusting the pd repulsion and 
combining calculations with and without the Ad electrons 
in the core of the In pseudopotential 0, • 

The OEPx calculations in the present work were 
performed with the plane-wave, pseudopotential code 
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par am. 
unit 
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CO 
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E 9 
eV 


Ax 

eV 


zb-GaN 
zb-InN 


4.50 
4.98 


3.07 
0.53 


wz-GaN 
wz-InN 


3.181 
3.533 


5.166 
5.693 


0.377 
0.379 


3.32 
0.72 


0.029 
0.067 





TABLE I: Lattice parameters, Go Wo band gap E g , crystal 
field splitting Ai for zinc blende (zb) and wurtzite (wz) GaN 
and InN. 



sfhingx [l8j. while for the Go Wo calculations we have 
employed the Go Wo space-time method [19j in the gwst 
implementation |2fl l2lj . Exact-exchange pseudopoten- 
tials [z| were used throughout and the cation d-electrons 
were included explicitly |TTL ITlj . The calculations were 
performed at the experimental lattice constants j2^| 
taken from Ref. [24| and reported in Tab. [I] Conver- 
gence to within 0.05 eV in the quasiparticle energies was 
achieved for a 4x4x4 (4x4x2) k-point sampling in the 
zinc blende (wurtzite) phase and a plane-wave cutoff of 
75 (65) Ry for GaN (InN). Unoccupied bands up to 55 
(45) Ry were included in the calculation of the polaris- 
ability for GaN (InN) in the OEPx as well as in the G W 
calculations. 

The Go Wo band gap E g and the crystal field splitting 
Ai for zinc blende and wurtzite GaN and InN are shown 
in Tab.HJ Since Go Wo falls into the realm of perturbation 
theory it is not a priori clear if the quasiparticle correc- 
tions are positive or negative. While LDA based Go Wo 
calculations generally open the band gap from the un- 
derestimated LDA value we observe here that the OEPx 
gap of 1.0 eV for wz-InN closes to 0.72 eV after applying 
the G0W0 corrections. This is not unexpected since the 
dielectric screening inherent to Go Wo but not to OEPx 
counteracts the exchange effects. 

Figure n shows the Go Wo bandstructure for wurtzite 
InN and the conduction band of wurtzite InN and GaN 
(circles). To make contact with experimental results 
we use an analytic expression for the conduction band 
around the T-point 



E c (k) = 
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2m Q 
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derived from a four band k • p model, neglecting spin- 
orbit splitting. Here mo is the free electron mass and 
E^ = E g ± Ai. The parameter E p is related to the 
optical matrix elements between conduction and valence 
bands. Since it is the only unknown in Eq. ^ it has 
been determined by fitting to the Go Wo conduction band. 
The small reciprocal lattice vector spacing required for 
an accurate fit is easily realised in the GoW) s pac e-time 
method by means of Fourier interpolation |20j | . 

Using the values of Table[]]fdr Ai and E g we find that 
the quasiparticle conduction band of InN is well described 

by Eq.Hand £?J = E£-=9.0 eV (red solid line in Fig. [TJ. 
For GaN the transition matrix elements differ for differ- 
ent directions and we obtain Ep =14.9 eV along the c-axis 
and £^- = 13.6 eV in the basal plane (black solid line). 
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FIG. 2: For wz-InN the effective mass as a function of carrier 
concentration, deduced from the Go Wo calculations by means 
of Eq. 0and|5| (red line), agrees well with experimental mea- 
surements (symbols) and the k-p fit of Wu et al. (black line) 
based on the experimental data |25ll26| . (Figure adapted from 

it) 
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FIG. 3: The Burnstein-Moss effect deduced from the Go Wo 
band structure (red line) for wz-InN reproduces the experi- 
mental trend (symbols) very well and is also close to the k-p 
curve of Wu et al (black line) |25l |26| . Assuming parabolic 
bands (dashed line) overestimates the Burnstein-Moss shift. 
(Fig ure adapted from Ref. the following symbols differ: • 
123 |23 A 0) 



While the conduction band is well described in the ef- 
fective mass approximation for GaN (blue dashed lines in 
Fig. ^) it exhibits a deviation from the parabolic shape in 
InN. Wu et al. have recently reported a dependence of the 
conduction band effective mass on the free carrier con- 
centration in InN pM |26| , which indicates a pronounced 
non-parabolicity of the conduction band. In the spherical 
band approximation the momentum effective mass 



m 



toq dE c (k) 



h 2 kp dk 



(2) 



k—kj. 



can be translated into a carrier concentration depen- 
dent effective mass using the free electron relation kp = 
(di^n) 1 / 2, , where n is the density of electrons in the con- 
duction band |2|| ■ Figure shows the effective mass 
of wurtzite InN as a function of the free carrier concentra- 
tion. The ab initio prediction (red line) extracted from 
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TABLE II: Conduction band effective masses m* , Luttinger (71-73) and A parameters as well as values for E v (zinc blende) 
and Ep and E p (wurtzite) [29^ for GaN and InN obtained by fitting a 4 x 4 k-p Hamiltonian to the Go Wo band structure. 



our Go Wo band structure by means of Eq.Qand El re- 
produces the experimental results very well and closely 
matches the curve obtained by Wu et al. with an ex- 
perimentally deduced value of £^=10 eV [251 EbT] (black 
line) . 

Equation ^ can also be used to calculate the shift of 
the direct optical transitions (E g (n) — E c (n)—E v (n)) to- 
wards higher energies in absorption measurements upon 
increasing electron concentration in the conduction band 
- the so called Burnstein-Moss effect. Contributions from 
the electron-ion and electron-electron repulsion at high 
electron concentrations are accounted for following Wu et 
al. [2^, |2|j . The Burnstein-Moss shift calculated in this 
way fits a wide range of experimentally reported measure- 
ments very well, as shown in Fig. [31 and agrees well with 
the curve deduced by Wu et al. from their experimentally 
determined values of E g and E p . Neglecting the non- 
parabolicity of the conduction band (blue dashed line) 
worsens the agreement with the experimental results. 

For device simulations k-p models have been well es- 
tablished. However, so far most of the parameters enter- 
ing the k-p Hamiltonian [30ll3l| for GaN and InN have 
not been conclusively determined by experiment Q. We 



have therefore used our quasiparticle band structure to 
derive a complete set of Luttinger and /^parameters for 
zinc blende and wurtzite InN and GaN [33] ■ The param- 
eters [2{j are shown in Tab. [nj 

In conclusion, we have carried out Go Wo band struc- 
ture calculations based on exact-exchange density func- 
tional theory for InN and GaN. From these results we 
have derived key electronic quantities (band gaps, carrier 
dependence of the effective mass, k-p parameters). For 
InN, a band gap of w 0.7 eV is found supporting recent 
experimental observations. For all investigated parame- 
ters an excellent agreement with experimental data has 
been observed, indicating that the wide interval of exper- 
imentally observed band gaps can be largely explained by 
the Burnstein-Moss effect. 
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